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ABSTRACT 

The Economic Order Quantity (EOQ) model is extended to 
a deterministic, two-echelon model with an arbitrary number 
of activities on the lower echelon. Two variations of the 
model are developed using minimization of time-average cost 
as an objective. A no-stockouts-allowed case is examined, 
and a method for finding the optimal solution is developed. 
A backorders-allowed model is derived and partially solved 
here in general. A full solution is presented for a re- 
stricted range of lower echelon parameter values. Examples 
of the no-stockouts-allowed model are given and solved. 

The solutions from this model are compared to those derived 
assuming the activities operate wholly independently. 
Significant potential reduction in variable time-average 
cost through the use of this model is demonstrated. 
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I. INTRODUCTION 



A field of primary interest in Operations Research, 
especially as applied to supply support problems, is the 
study of inventory models. Inventory theory is the appli- 
cation of mathematical optimization techniques to the problem 
of deciding when and in what quantities certain goods should 
be purchased to meet the demand of customers while minimizing 
time-average cost [1]. 

The fundamental inventory model, usually called the 
economic order quantity (EOQ) or Wilson model, presumes that 
all costs and demand are fixed over time, that resupply 
occurs after some deterministic leadtime following placement 
of an order, and that demand will be met without incurring 
shortages of stock, all of this for a single activity hold- 
ing only one good or commodity (see the sketch below) . 





PurcV^ase Q. 


Mi 






CoMst<vv^-t 
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More complex and realistic models have been developed by 
relaxing some or all of these assumptions [1]. Nevertheless, 
the majority of the problems that have been solved concern 
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finding optimal policies for single activities since multi- 
echelon problems are more difficult to formulate and solve. 

For many practical applications inventory models for a 
single activity suffice. But for some large retail firms 
and military supply systems which have centralized procure- 
ment for a system of retail outlets or stock-points, a 
mathematical model which accounts for the multi-activity, 
multi-echelon nature of the inventory system is more appro- 
priate. A classic principle of optimization theory is that 
any allocation of resources for a multi-component system 
based on optimal solutions for each component operating 
independently can be no better than the solution found for 
the system as a whole by subordinating the roles of the 
independent components. It is the author's belief that most 
mult i- echelon inventory systems tend to operate sub- op timally , 
because the above-mentioned independent component approach 
is used to solve for stocking policies. 

This paper presents a simple, deterministic, multi- 
echelon model which is optimized on a "system" basis. The 
model is an extension of the EOQ with two echelons, as 
depicted below. 
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The upper echelon consists of a single activity which could 
represent a centralized procurement and warehousing facility. 
The lower echelon has m independent retail outlets or stock 
points which replenish from the upper echelon and meet the 
demands of the customers outside the system. A minimum 
time-average cost solution for the system, presuming no 
stockouts are allowed and subject to the system cycle length 
being an integer multiple of the cycle length of each lower 
echelon activity (possibly a different multiplier for each 
activity), is derived in Chapter II. In Chapter III, an 
algorithm is developed, to find the best choice of integer 
multipliers. The case in which customers' demands may be 
backordered is taken up in Chapter V. Some examples of the 
no-stockouts-allo\>fed model solution are presented in Chapter 
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IV with an analysis of the cost benefit available in the 
use of this model over the use of independent EOQ models 
for the activities. 

The work in this paper is basically an extension of 
the two-echelon model with one activity at each echelon 
sketched in [2] and presented in greater depth with the 
finite production rate and backorders cases in [3] . 

Several authors have attacked special cases of the multi- 
echelon problem. A probabilistic model for repairables is 
presented in [4]. A probabilistic two-echelon model is 
solved in [5] using a dynamic programming technique. 

The deterministic models in [2] and [3] , while they are 
fine in themselves, solve problems of completely specified 
size. The stochastic models of [3] and [4] are not multi- 
echelon extensions of the probability models presented in 
Chapter Four of [1] . This paper is a natural two-echelon 
extension of inventory theory's fundamental model with an 
arbitrary number of activities at the lower echelon. It 
is hoped that it also provides a starting point for further 
development of more complex deterministic models and more 
general stochastic models such as extensions of those in 
Chapter Four of [1] . 

The author presumes in his presentation that the reader 
is familiar with the EOQ model and its backorders- allowed 
extension. 
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II. THE NO-STOCKOUTS -ALLOWED MODEL 



A. MODEL FORMULATION 

Suppose that there is an inventory system consisting o£ 
two echelons of inventory, performing the procurement, the 
holding, and distribution of a single good. The upper echelon 
is a single activity which orders and buys the good from an 
outside source and holds it in inventory for further distri- 
bution within the system. The lower echelon has m independent 
activities which order and resupply from the upper echelon 
and hold inventory to meet the demands of their customers. 

In setting up the cost equation for the system the following 
are assumed: 

1) All orders are filled immediately (a deterministic 
lead time has no effect on the optimal policy) ; 

2) all demands are met immediately; 

3) the demand rate is constant and continuous over time 
but may be different for each lower echelon activity; 

4) the goods are purchased only by the upper echelon; 
and the price is constant over time; 

5) ordering costs are constant over time; 

6) inventory holding costs are products of the annual 
interest rate and the purchase price of the good 
(following [1]); and 

7) the cycle length (time between orders) will be constant 
for each activity; and system cycle length will be a 
positive integer multiple of each lower echelon activ- 
ity's cycle length; the number of complete sub-cycles 
for each lower echelon activity may be different. 
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The sketch below is a schematic of the material flow in 



this inventory system. The construction of the mathematical 
model of the system parallels that of the EOQ model. 




Throughout this paper whenever a summation sign ( T ) is 
used, it is assumed that sum is taken for the index variable 
ranging from one to m unless it is specifically stated 
othen\fise . 

Activity 0 is the upper echelon activity. Activities one 
through m are the lower echelon activities. The subscript of 
a variable identifies the parameter or variable with its 
associated activity. The A^, and are, respectively, 

the ordering cost, interest rate, and consumer demand rate 
for the i^ activity. These, along with C, the unit price 
of the good, are the system parameters. The Q^, T^, and n^^ 
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are, respectively, the ordering quantity, the cycle length, 

1" h 

and the number of complete cycles over for the i^ activity; 
Q. and T. are continuous variables while n. is a positive 
integer variable. 

Let H. be the total costs of running the inventory system 

for one system-cycle, T^. Since each cycle is comprised of 

th 

n^ sub-cycles for the i^ activity, total ordering costs are 
. Total purchasing cost is merely CQ^, for the 
goods are bought only at the upper echelon. Let Bj^(t) be 
the quantity of inventory on hand at time t. Then the total 

»v> 1 

holding costs are ^ IlC J BiCt)clt . Putting the parts 
together ; 

To 

Z = Ao ^ + £ liCj BiWjt. 

>— O i. O ^ 

A typical realization of the on-hand inventories within the 
system over one cycle is shown in the following sketch. 
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The assumptions of the model result in: 



T = n. T. for all i; 
oil 

for all i; and 



• I nj Qi . 

Now consider liC J BiW Jt 



( 2 - 1 ) 



for Since 



the n^ sub-cycles are identical and the area under the on- 

hand curve for one sub-cycle is just that under a triangle, 

T: 

QiT, 



= nJ-,C | Bi(.t)clt= /ivIiC 



For the upper echelon activity Bq (t) can be decomposed 
into the sum of the inventories held for each of the m lower 
echelon activities. If (t) is the inventory held for the 
i^^ echelon, then B^ (t) = 2 B^ (t) where 



{n-1 - Int[^] 



Q; 



and is 



the greatest integer less than or equal to X. 
Then 









jco 

Summing over i gives 

To 



I Bo tt'id t = I^Z . 



Z- 



It follows that the total costs per cycle are 

A.- + C(5. + 21 i r.cGii 4 z, i r.ca.Aa-it-i)T^. 
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B. MINIMIZATION WITH RESPECT TO Q 

To compute the time-average cost for the system, divide 
the cost per cycle by the cycle length, T^ = T^. Thus, 

K, ^ lirCQ; V (i-Z) 



Although this is a function in 2m + 1 continuous vari- 
ables, it can be reduced to a function in one continuous 
variable by using the equations (2-1) that arose from the 
model's assumptions. Since Q. = D- T- and T„ = n- T. = 

^ ^111 oil 

n^ Tj for i,j = l,...,m; equation (2-2) becomes 



K-- 



A,D, y ALyiL-D. 
0.(5, * ^ 



+ c + 



T IlCDlv^.Q^ y T.C».g(nL-l) D: 
^ ar>, ni ^ 2D,ni 



(2-3) 



In this form, K will be referred to as K(Q, , n, , . . . , n^^^) . 

For a given choice of n, , . . . ,n K(Q,) will represent equation 
(2-3). 



Taking the derivative of K(Q,) gives 

iJ< - . AeD, ^ y A.,AvP. j. TyC D^n , j loC r\.Crty-l)D: 

dQ. ry,Q} ' ^ ^ 2D,ni ^ ^ 2D, r?y ^ (2-4) 



and 



Now , 






O implies 




' 2 {Ao^KA^J 


I 






L ‘ 1 J 







24.D, 
n. Q.3 







o 



(1-5-) 



since all of the parameters and variables are positive. 

A 

Therefore, Q, gives the minimum of K(Q,). 

Substituting Q, given by equation (2-5) into equation 
(2-3) gives 
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> % 



L'l'fo) 



Let (Q, , rt, ,-.. , ) be an optimal value 

of (Q,jn., , ) . Then 

Theorem 1-1 : K(Q> »''>) '"i " K (Q< » ’^' > ** * j 

Proof : Since K(Q<^ j 

K KCQ., . Now, for each fixed 



- K for all Q^. 

Then K £ K (Q«>nr)— >«« ) for all 

Thus, K C^' Finally , 

KCQT,nr,...,nJ^) < ^ K (Q^Ar.-.-.O) which 

implies the result. 

y\ 

Corollary 1-2 : = <3i • • • > ^<^ ) . 

C. CLOSING REMARKS 

"h 'k k 

Given Qi ,n the remaining components of the 

optimal operating doctrine can be derived. Let starred 
variables represent the optimal values of the variables. 
Then 



rr -- 




• 

) 


(2-7) 


Ti'- 


A, -|- * 
A; ' • 


*p6r c^n L 


(2-8) 


M 




•for all 1 J 


(2-9) 


r/= 


nrv 


0. ni 


(2-10) 
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Q.' = I VI* Qi* . 



( 2 - 11 ) 



In this chapter, the no-stockouts -allowed model has been 
developed, and the values of the continuous variables which 
minimize the time- average cost function (given any set 
, . . . , n,v, of the integer variables) have been found. 

In the next chapter, a method for finding the optimal values 
of the integer variables is derived which will, by using 
equations (2-7) through (2-11), give the optimal policy for 
operating the inventory system. 
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III. FINDING OPTIMAL n IN THE 

NO- STOCKOUTS- allowed MODEL 



Chapter II developed the no-stockouts-allowed model and 
gave in equation (2-6) the minimum value of K for a choice 
of n, ,...,nni. In this chapter a method for finding optimal 
n|,...,n,y, is given and, for the case A^ = A, for all i, an 
optimal search algorithm is stated and proven. 



A. THE GENERAL CASE 

Let n be the m-tuple (n,,...,nj^) of positive integers. 
Let K(n) be the optimal K for each choice of n as given by 
equation (2-6). Let n* be the optimal value of n. 

Minimizing K(n) could involve considering a countably 
infinite number of alternatives for n*. Another potential 
problem in seeking n* is that there may exist a sequence 
of n's, call it jn .1 , such that the Thus, 

— C — J J j-»oo 

even though K(n) is positive and, hence, bounded below, it 

may not have a minimum.^ 

Fortunately, neither of the above problems arise with 
K(n) . Theorem 3-1 and its corollaiy . given below show that 
n* lies in a finite, one-sided neighborhood G of the point 
(!,...,!) whose boundary can be computed easily. Hence, n* 
can be found by considering a finite number of candidate 
points n. 



As an example of this peculiarity, consider y(n) = 
for n = 1,2,3,... Clearly, y(n) > 0 for all n; yet 
Hence, there does not exist a finite n* such that y(n*) ^ y(n) 
for all n. 
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Theorem 3 



- 1 : For each j 5 , there exists an Nj 



such that K(n) > whenever n^ ^ regardless of 

the choice of n^, i j . 

Proof : Let 

U3>= 51 I KCQ) - ClDif 

Note that LCh^> L(i, •••,!) only if 

Let be arbitrary 

but fixed. Let 



0< = Ao 

c4j 

0^ I ^ X I D L 



and 



^''' h. 

Then o( > O and ^, > O . Then 

L<».. 

: Cn^-OA.I.Bi * A,t,D, 

+ |c</3, *o( vc< 



Let 



Since LCO ) - («i) > 0 , L CQ ) > C ^ . 

b - Ul, ->i) , and let 

Ni- 



Let 



where Int [•] is the greatest integer function, 
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Since MCn^) is a strictly increasing function in 

N-, M(n.) > M(N.) for n- > N . . But 
3 3 J J j 



N.> 



b — Aj r j Dj 

Ajio D; 




and 



MCMO > ) A-.nDi A-.XiDi 

= b = LCt,...,t). 



Hence, for all n^ :> and regardless of the value of 
the n^ , i # j > LCQ')> 

Then for the same conditions on n, K(n) > K(l,...,l); 
and the proof is complete. 

Theorem 2-1 shows that the neighborhood 
Gc - [ - ’ i - ^ is a hyper box in 

the positive integer lattice in m-space with one corner at 
the point (1,...,!). The following sketch shows what G 
looks like if n = (n i, n^,) . 
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An inunediate implication of theorem 3-1 is the following 
corollary that n* lies in the bounded hyper box G. 

Corollary 3-2 : For K given by equation (2-6) 

t. . 

Proof : Either n* = (1,...,1), or it is not. If 

n* = (1,...,1), then n* = G. Assume now that 
n* ^ (1,...,1). Then K(l,...,l) > K(n*) . Now 
for all n ^ G, 

K(n) > K(l,. . . ,1) > K(n*) 

which implies that n* is not an element of the complement 
. of G. Hence, n* £ G . In either case n* £ G. This 
completes the proof. 

The previous theorem and corollary show that in the 

general case n* and, hence, the optimal policy ( (2t , Ti. ) 

can be found in a finite number of steps by an exhaustive 

search. The exact number of steps required is equal to the 

number of positive integer lattice points in G; namely, 

TT For example, suppose that m = 5 and 

1 = ' 

= 11 for i = 1,...,5. Then the ntimber of iterations 
required to find n* is 100,000! 

B. THE CASE Ai = FOR ALL i 

The balance of this chapter is devoted to developing an 
optimal search algorithm to find n* for the special case of 
having the same ordering cost at each of the lower echelon 
activities; that is, A^ = for all i. The method of search 
• is demonstrated by the following example in which n = (n , , nj^) . 
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Suppose N, =5 and N 2 , = 3. The search begins by com- 
puting K(l,l) [ actually L(l,l)j and saving (1,1) as the 
first candidate for the optimal solution. The two adjacent 
lattice points in G, (1,2) and (2,1) are compared to find 
out which one gives the smaller value for K. Say that 
K(l,2) >K(2,1). Now K(2,l) is compared to K(l,l) to find 
out whether (2,1) will become the updated candidate for n*. 

If K(2,l) < K(l,l), then (2,1) becomes the updated candidate. 
Whether the candidate changes or not, the search continues 
from the point (2,1). Now the two points adjacent to (2,1) 
for which either n, or increases (that is, (3,1) and 
(2,2)) are compared in the same manner as were (1,2) and 
(2,1). Say that K(2,2)< K(3,l). Then K(2,2) is compared 
to the candidate for n* as before, and either an update 
occurs or it does not. In either case the search continues 
from the point (2,2). The search terminates when either 
n, = N, or = N^, and, at that point, the coordinates of 
the last candidate give n*. The diagram below shows how 
the above search began and how it may have continued. 




on Path 


CAudtdafe 0 


c.o 


Cl, 1) 


u»o 


(•>«) 


(i.i) 


C2.,^) 


(3,2) 


C2.2-) 


(3,3) S+op. 
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1 

I 



I 



The algorithm, itself, will be stated at the end of 
the chapter. Some theoretical groundwork is required to 



prove that this method of search works. 

In the following development, let g (n) be a positive 
increasing function in each n^ defined over the positive 



integer lattice in Euclidean m-space. Let n^ = (1, 



, 1 ). 



Let element in a sequence of n * s 

which track the progress of the search through the positive 
integer lattice. Then let 

and let n^j^ be such that 



Then the sequence [ S, j | traces the coordinate-wise, unit 
step-sized path of steepest ascent in g(n). 

Define h Ql) by 



Thus, h is defined on all positive integers and corresponds 
to the value of g (n) along the path of steepest ascent. 



Now, when = A j for all i, K contains a function 
g(n) for which it is important to find the coordinate-wise 
path of steepest ascent. Without justifying its existence 
at this point, the particular g(n) will be defined and a 
property of its associated function h will be noted. 
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Theorem 5-5 : Let g(n) = ZL where gi(0) = 0 

and gi(n^) = ^ Suppose that 

I > 0 andX> > 0 for all i. Ifh(J^) is defined for g(n jj ) 
as above, then h (JJ ) g (n) for all n such that 



Proof : 



Let 




where 



t "*1 

Let h (0) = 0. Finding the j( — element of the sequence 
is equivalent to picking the ^ largest elements 
from the set and letting n^j^ equal one plus the number 

of elements picked from • Then h ( ) = g (n^^ ). 

For any other n such that s ( L - yvi. , 3C^'i is 

found by picking the largest n^^ - 1 elements , taking 

their sum, and then summing over i. Since the sum of the 
largest elements in must be at least as large as 

A ..-I 

I. ZT (j — largest element in ) such that 

0 - ' 

r (l - m; 

and the proof is complete. 

Consider l-» (n) with = A for all i. 

[A.+ A,L»';j[l.£D;<- I. 

Let - Iq and suppose for the present that > 0 

for all i. Let the set J = . Then 

LC0)= [Ao+ A,L n_|[roZ.t>i + XivD;.- } . 

The last expression in the second set of brackets is g (n) 

in the previous theorem. 
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The following theorem completes the proof that 
coordinate-wise steepest descent works for L (n) and, 
hence , for K ( n ) . 



Theorem 5-4 : Let n^ be 



min 

n 



L (n) = 



V] tW 



Let}) . 



defined as before and let 

T -C niiTi j \ ^ 

. If n L ( n ) exists , 



then 



Proof : First note that theorem 3-1 and its corollary 

showed that HjLO) exists. Let a ^ = Aq , a^ = A, , and 
aj = LTv.B; . Then 

XiClO^ - [a , + I etj - 

Clearly, 5 wm’a * Suppose 

"* !S t'w 

I_j (.»' ) < HcH'i . Since the minimum exists, there 

« « tW 

is an n i W such that jL(vi) <»^i/).XjCO) . Let S ■= . 

— i Vv ' 

A 

Then, from theorem 3-3, - 3C0^) * 

It follows that 

L (0 ) = { <* . * 4- 2 ■’;} - 3 wij. [ 4 ,* (is - 3t2 >1 



This contradiction gives the result. 

Corollary 5-5 : ^ KCO^ = *<(.!? 

Proof : See the proof of theorem 3-1. 

Theorem 3-3 assumed that > 0. The following theorem 
shows why this assumption is not restrictive. 
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1 . 



Theorem 



5-6 : If Ij ^ 0 for some j, then nj 



* = 



Proof : Suppose that nj* > 1. Let n* be such that 

XiCO*) =■ • Now 

■{ )( T XvDl XiDil 



Let n - (n nj 1 , . . . , ^ ^ Since n . < n , 

J J ’ 

(a , + aj^l, n^) < (a , + aj^ n" ) . Since I-j ^ O , 

.< and (I„ r D, 

Then L (^) < L ( n*) . This contradiction implies n j * = 1, 
and the proof is complete. 

For the special case of for all i the preceding 

theorems have proven that the following algorithm will 
provide the optimal solution. 



Algorithm 

(1) Let J = I'- ' . If J = ^, stop; n^* = 1 for 

all i . Otherwise , set n = (1,...,1). Compute L* = L ( n ) 



and N. = Int [ + i] + i. get n = 

1 At Xo D*^, “ 

and go to step (2) . 



( 2 ) 



Let j £ J be that index such that 









for all i £ J. Increase n. by one. If n. = N . , stop; 

3 3 3 

n* = n. Otherwise go to step (3) . 

A 

(3) If L (n) < L* , set L* = L (n) and n = n. Otherwise, 
make no changes. In either case, return to step (2). 
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After n* is found, then Q^*, and K* can be computed 

as shown in Chapter II (see equations (2-6) through (2-11)). 

Earlier in this chapter it was shown that for an exhaustive 
search over G for which m = 5 and = 11 for i - 1,...,5, the 
number of iterations required would be 100,000. If this 
algorithm could be employed (that is, if = A^ for 
i = 1,...,5) for this case, then the maximum number of 
iterations required would be 50, an impressive reduction. 
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IV. EXAMPLES OF THE NO-STOCKOUTS-ALLOWED MODEL 



This chapter is dedicated to some examples which, it is 
hoped, will enhance the understanding of the model and will 
aid the reader in following the use of the algorithm 
developed in Chapter III. 

A. AN EXAMPLE WITH LOWER ECHELON SYMMETRY 

A two-echelon, three- activity inventory system is 
operating with the following parameters: 

Aq = $500 ; Ai = Aj^ = $1 ; 

Iq = 0.1/year; ^1 “ ^2, “ 0.05/year; 

C = $10/unit; D = Dj. = 1,000 units/year. 

The aim is to find the optimal operating doctrine. The 
following sketch illustrates the system. 
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The algorithm in Chapter III terminates at step (1) 
since = I 2 ,. Equation (2-10) gives 

K*(l,l) = $21,002 (or $1,002 in variable costs 
alone), and equation (2-5) results in 

Q * = 1,002 units. 

Equations (2-11), (2-12), (2-13), and (2-14) yield the 
remainder of the optimal operating policy: 

T * = 1.002 years; 

T*=T* =T** 

^o ’ 

Qi* = Qi*; and 

Q * = 2,004 units. 

^o 

Had the inventory policy been derived using the standard 
EOQ model for each activity, the following would have resulted. 
Since activities one and two are identical, 

= 63.24 units; and 
Ti* = = .063 years. 

This results in a time-average variable cost for each of 
the lower echelon activities of 

Ki* = $31. 62/year. 

To find the time-average cost for the upper echelon recall 

that T^* = T 2 * and, hence, that the upper echelon must 

supply [ Qi* Q<i*j every T^ * units of time. The holding 

and ordering costs divided by T^ for the upper echelon are: 

Ao^ ^ loC 

Z ~ K(«i). (4-1) 
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(4-2) 



Since Q^* = Q^*, equation (4-1) becomes 

K= - r»CCn,-0Q* . 

Treating nj_ momentarily as a continuous variable, it is easy 
to see that equation (4-2) is convex in n^^. Minimizing with 
respect to n^ gives two candidate integer solutions for the 
minimization problem, n^^ = 11 and n^^ = 12. Of these, 
nj^ = 11 gives the lower cost value for the problem parame- 
ters. Then the total variable time-average cost in this 
second analysis of the inventory system is 

+ r.CCn'-l) <2* * 1Y.I 

= $1,413.40. 

Comparing the second analysis to the first, the second 
policy yields a cost which is 41% higher than that given by 
the optimal policy for the model of Chapter II. This seems 
only natural since the model of Chapter II allows for 
cooperation among the various activities of the inventory 
system. 

B. AN EXAMPLE WITH LOWER ECHELON ASYMMETRY 

The following example has some differences in the 
parameters of the lower echelon activities and the resulting 
cost difference between models is more dramatic than that 
in the previous example. 

Consider a two-echelon, three-activity inventory system, 
as in the previous example, with the following parameters: 
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Aj^ = $100; 



= $500; 



A, = $1; 



= 0.1/year; “ ^2 “ 0.05/year; 

C - $10/unit; = 100 units/year; 

and = 1000 units/year. 

Again, since 4 0; n,^* = n^^* = 1. In this case 

the total variable time-average cost is 

K* = $813.08. 



The optimal operating policy is: 

Qj* = 147.83 units; 

= 1,478.3 units; 

T. * = T. * = T * = 1.478 years; and 

Q * = 1626.13 units. 

^o 

As in the first example, a policy can be derived assuming 
independently operating EOQ activities and the variable costs 
can be compared. 

The EOQ operating policies for the lower echelon activ- 
ities are: 

Q^* = 200 units; 

Qj,* s 63.24 units; 

T^ * = 2 years ; and 
T 2 * = .06324 years. 

This yields time-average costs for the two activities of 
K^* = $100/year , and 
= $31. 62/year. 
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There are two alternative policies for the upper echelon to 
employ. One is to order for both lower echelon activities 
on the same cycle. This leads to the constraint 
n^T^ * = n<^Tj^*. The smallest values of n^_ and n^ which solve 
this equation are n s 1,581 and n<^ = 50,000, and they 
result in a policy with ridiculously high costs for the 
upper echelon. The other policy is to order for the two 
activities on independent cycles, and it is this second 
policy that is derived. The variable cost for the upper 
echelon as a function of time is 

= AoT*'t[^;T7] + A9lntLixTj‘l + holding costs (t) . (4-3) 
Dividing equation (4-3) by t gives 



Kw- 



Ao I t J.* ] AoTnt 



+ ave. cost^Ct) . ^4-4) 



As t grows large, Int *] = X.(t) , a large integer; and 

n ^ 1 ^ ' '1 

t " "i " 



Similarly, 



t 



L 



Average holding costs, by an argument similar to the one 
above, stabilize at 

i.c (n.-Da* ^ T^c ci.-iLQ» ^ 



Then |im Kit)- K , where 
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K 







Ac 



(4-5) 



Ai, + T.c TeCK-DQ^: 

A — — 



= K(n^,nj) 



Equation (4-5) is convex in n and r\. 2 ^ if n^ and nj_ are 
treated as continuous variables. Minimizing K with respect 
to n^ and nj^ gives n^ = 1.58 and n^^ = 15.8. Four alternative 
solutions for integer values of (nj^,n^) are (1,15), (2,15), 
(1,16), and (2,16). Of these four, the choice that offers 
the smallest value of K is n = 2 and n 2 ^ = 16. For the 
upper echelon alone, then. 



K(2 , 16) = $1,193.45. 

Adding to this the time-average cost for the lower echelon 
activities gives a system time-average variable operating 
cost of $1,325.07. 

The second operating policy yields a system time-average 
cost which is 63% higher than that for the optimal policy 
for the model presented in Chapter II. With asymmetry in 
the problem, allowing the activities to cooperate pays even 
greater dividends. 



C. AN EXAMPLE WHICH EMPLOYS THE ALGORITHM IN CHAPTER III 
This final example shows more details of how the 
algorithm of Chapter III works. 

Consider a two-echelon, three-activity inventory system, 
as before, with the following parameters: 
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A = $5 ; 
o 

= 0.05/year; 
C = $10/unit; 



= $5; 

^ 1 “ ^ 2, “ ® • 1/year ’> 

Di = 100 units/year; and 
D^, = 1000 units/year. 



For = 
N ^ = Int 

N 2 , = Int 
Since for i = 



n^ = 1, L(n) = $1650. 
1650-50 



[ 
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+ 1 ] + 1 = 66 , 



rl650“500 j. -1 "I ^1 _ £ 

t^50 + 1] + 1 - 6. 



1 , 2 , 



=5/2 and 



UDg - 



50/2 ; 



Then 

and 



set = 2. Now n = (1,2) and L(n) = $1700. 

Then n remains (1,1) and L* remains $1650. 

The table below shows the progress of the algorithm. 

In this example nj_* = nj^* = 1. The optimal policy could 
then be computed using the equations at the end of Chapter II. 
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LCO) 


L* 


A 

n 


1 


i 


1 (o So 


1 G 5 0 


a.i) 


i 


X 


iTOo 


1 C,5o 


(1.1) 


1 


3 


I'l 


I&5 0 




1 


4- 


aits’ 


(450 


Cl.i) 




4 


^450 


(45 0 


d.i) 


X 


5 




/450 


(1,1) 


3l 


ip 


Stop, Sl^lce . 
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V. THE BACKORDERS -ALLOWED MODEL 



A. MODEL FORMULATION 

This chapter extends the model developed in Chapter H 
by allowing the lower echelons to run out of stock and to 
backorder the shortage until the next order is received. 
For each unit backordered a penalty cost is assumed to be 
suffered. 

The sketch below shows the net inventory during one 
cycle for one lower echelon activity. For activities one 
through m a new variable, s^, the backordered quantity for 
one cycle, and a new parameter it^, the shortage penalty 
cost, are introduced. 




(5-1) 



The total costs for one cycle are given by (5-2) 




+ Z (holding cost for lower echelon activity i) 
+ Z (shortage penalty cost for lower echelon 



activity i) . 



(5-2) 
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The holding cost is the product of I^C, the unit holding 
cost per unit time, and the area under the on-hand inventory 
curve over n^ periods. That is: 



The shortage penalty cost is the product of the unit 
penalty cost, and the total number of backorders incurred 
over n^ periods. That is, n^ir^s^. Thus, equation (5-2) 
becomes : 



"2 - Ao L^‘Ai + C + loC L"^v 






a-. 






(5-3) 



Dividing equation (5-3) by T^ gives the time-average 
cost for operating the inventory system. 



K 



n. Q. 



C IDi 



+ T.C L 



Cy>..-0>^.'D;Qi 
5 y)'^. 



z 



A V. 

(a, 



+ 




o.DiQi 

2nvD. 



S -, -H 



Sv’'lO tl>i I j. T 

u:e:s: s ^ 



Try’S, n-.D, 
n.Q I 



(5-4) 



Equation (5-4) will be denoted by K(Q^ ,s^ , . . . ,Sjj^,nj_ , . . . ,n^^) . 

B. OPTIMAL OPERATING DOCTRINE FOR GIVEN n ,...,n . 

m 

A full statement of the problem faced in minimizing 
time-average system cost would be: 

minimize K 5,)... jSm, O, j... ) 

subject to Qo- X'^lQl Ql-DlT, -fcrr l- l,--> > 

1 aT^<i 
O -foTT 

and n.X - n-T, = To -for 
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Equation (5-4) presents K with the equality constraints 
already satisfied. Full justification for this simplifica- 
tion is not presented here, but the reader is invited to 
consider what information the Kuhn-Tucker conditions would 
give in the case of the equality constraints and how this 
would then be used to simplify the other first order condi- 
tions. Consider only the problem posed in minimizing K 
given by equation (5-4). 

minimize K(Q^ ,s^ , . . . , . . . jn^^) 

subj ect to Q 0 ; 

n^ 1 and integer. 

Disregarding the variables n^,...,n^ for the time being, 
notice that the constraint set is convex. If K itself is 
convex, then the entire problem is a convex programming 
problem with a global optimum. The issue of the convexity 
of K will be discussed later. 

/V A 

Let , Sj_ , . . . , Sjjj denote the values of Q > • • • >5^^^ 

which minimize K for a given choice of n^,...,n . Let 
Ql* , . . . , Sjjj* ,n^*, . . . ,n^^ be the optimal values for all 
the variables. By arguments similar to those given in the 
proof of theorem 2-1 the following can be shown to be true: 



Theorem 5-1 : K(Qr,srv..jSm,n,%...,"« ). 



Now Q^,Sj^ , . . . ,s^ can be found by the calculus 



^ A A 

K(Qi ,Sj^ , . . . ,s ) is a global minimum for given n^,...,n 



m 
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if the first partial derivatives of K with respect to 
Q^,s^,...,s^ evaluated at Q^,s^ , . . . ,Sj^,rij_, . . . are zero 
and the Hessian of K is positive definite (so K is convex) 
for the specified n<,...,n and for all feasible values of 
Q^,s^ , . . . There are several cases, however, in which 

this approach to solving the problem fails because either 
the partial derivatives cannot be set equal to zero simul- 

A 

taneously for some ^ 0 and some s^ such that 

A A 

0^ St i. Qi for each i, or the Hessian of K is not positive 
definite . 

Taking the first partial derivatives of equation (5-4) 
with respect to Q ^ and s^: 



3j< _ _ AoD . j. T r T I ^ iTt AiDi 

5q, “ ^ D, - -Qx — 

TvC n.D: _ J_ y TtCs;^n:D» j_ y 7t>s->n,l) , 
^ Ao-.D. ^ “ Q?- «. 



(5-5) 



; and 



I-.c s-.n-tD, ^ Ttv 

35: = ^ • (5-6) 



Setting equation (5-6) equal to zero and solving for s^ 
yields : 



D- (n.Q.tvC - TL-wn-tD.) 

It C D, 



(5-7) 



Substituting this result into equation (5-5) and setting 
(5-5) equal to zero gives a quadratic equation for Q 
in terms of the system parameters. 



Bi + Ba = where 



(5-8) 
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and 




/ y n, I>L I k C 




J ^ i^kb. 


■ - 'i 




i 


L Xn-kt, ^ 






*4. A w ^ I 



loC 



^ an;D, 



4 - 



L 

I 



niDiK-^D, 

gn.x-.c 

n,I>iTwC 

a»Ki>i 
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Depending on the values of B and B 2 the method for 
finding Q^jS^ , . . . ,Sj^ varies so that a case by case analysis 
will have to be made. Note that it is possible that B^^ can 
be negative, zero, or positive, depending on the 
its numerator, 

Ao D , y iOl A tT), y TL^3)> 

^ ^ ^ 'O I A-j 3k v\ , ^ C, • 

Note also that if n^ = 1 for all i then B^ = 0; 

Bj^ < 0. Six cases will be considered below: 



CASE 




Bi 


ONE 


= 0 


= 0 


TWO 


< 0 


> 0 


THREE 


= 0 


> 0 


FOUR 


< 0 


= 0 


FIVE 


< 0 


< 0 


SIX 


= 0 


< 0 



sign of 



(5-9) 

otherwise , 



In the analysis of the cases, the following property 
of K will be valuable. For a fixed, positive Qj,, K is 
positive definite everywhere for s^^... ,s^^. If some 

of the s^ = 0 and Qi> 0, then K is positive definite every- 
where in the remaining s- ? 0. To see that this is true 
consider . 
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_ IiC 

as;‘SSi ” 



for i = j ; and 
for i # j . 



= o 

Hence, the Hessian of K in any non-void subset of the s^'s 

is a diagonal square matrix with all of the diagonal terms 

positive. This implies its positive definiteness and "the 

resulting convexity in the s^'s. Additionally, when the 

3 K 

value of s^ that renders = 0 for each i is substituted 

into K given by equation (5-4), K is then convex in Q^. 

This method for finding the minimizing values of 
,s^ , . . . , and Sjjj is justified in Section D. 

CASE ONE : = 0 and = 0. 

In this case equation (5-8) cannot be used to solve 
for Qj^ since it yields no information. Let > 0 be 

arbitrarily large but fixed. Then s^ is given by equation 
(5-7), and the positive definiteness of K in s^,...,s^jj 

A A 

implies that s^,...,Sjj^ provide a global minimum for the 
given Qj^. Now consider K as a function of at the value 

A A 

Si , . . . ,s . Its value is 

* m 

K = C L Bt • 

A 

Since K is independent of Q^, could be any positive value. 
Once it is selected, s^ could be computed using equation 
(5-7). 

A 

A problem arises, however, if is chosen such that 
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since equation (5-7) would give < 0 for some i. The 
problem is resolved by considering the following two state- 
ments of the problem. 

A 

Minimize K (Qi.s, ) . 

j max 

given Q, < ^ 

Minimize K(Q^ ,s^ , . . . ,Sj^) , 

^ . WAX TLit>i j 

given Qt < ^ ; and 



s^ ^ 0 for all i. 

Since the second statement is a restriction of the first, 
the best answer to the second can be no better than the best 



answer to the first. This implies that K(Q^,s^ 



,s ) for 
’ m-' 



^ . vn a. X ^v.^‘ 

liUw XX 



is greater than or equal to ClDi+ 



and that there is no advantage to choosing Q| that small. 

^ _ A ^ y TC y t 

Thus, Q, can be picked as long as Q, ^ i<i<m "iTc > 

A 

and s^ is then computed using equation (5-7). 

CASE TWO : ® j. ^ 0 > 0. 

A 

In this case Q can be computed using equation (5-8); 



namely Qi 



■ f- tel' 



*M«.X 



nLltilii 



Provided that n, TiC ") 



will be non-negative for all i; and the problem will be 
^ Wax 

solved. If Qi. < rt,XyC 



, then s^ < 0 for some i 



would result. Thus, for Qi < 
of solution is invalid. 



, this method 
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w^ax 



t)i 

Suppose that 

A A 

for s, that minimize K results in 



and fixed. Solving 



r I exprtsSionC5-*^)j + c iDl + 2]TC.lI>i+ %|loC Z 



= ~st“ j expression I 

I 

/N m <x X 

Additionally, suppose that Qj^ < i<Uvn n^T^C • 



Then 1< (Qo < minimum K subject to 

^ /s^ A 

Qt = Qi • In fact, as decreases from towards Qi_ , 
K decreases. Let 



rvia X 



niTCvDi 



^1 " ^ 

/V 

and let decrease from Qj^ to Q ^ = e^. Minimum K is taken 
on where s^ = 0 for at least one value of i. Let 



Yl 

Y., 



-[i ; Si = oj i 
= ^ i : Si > 0 



and 



Now as is decreased, the minimization of K given < e^ 
becomes 



minimize K 
Sc-- L* Yz 



subject to fixed; and 

Sj^ = 0 for it Yj^ . 

The functional form of minimum K is: 

\y — ■“ ^ ^ 1 , y n'c A c ^ I ^ iVcbiTt-i ) 

~ I n, ^ n, " In.IcC [ 






+ 



2L TL'cl>l -t 




I 



(n-v-l)n.I)L 

2n-,p, 






T niPL ) 



(5-10) 



The first bracketed expression remains positive and increases 
as the number of elements in Y^. decreases. The second 
bracketed expression also remains positive and increases as 
the number of elements in Y, increases. 
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Equation (5-10) remains a valid expression for minimum 
K in sj^,...,s^ given Qi fixed until is decreased to the 
point where equation (5-7) implies that another s^ = 0 for 
some i. Then the functional form of minimum K changes, 
because the sets Y, and change. Once Y^_ and Y^, have 
been modified to restrict any of the s^'s from turning nega- 
tive, equation (5-10) is again valid. This method of 
recursively formulating K remains valid until has been 
decreased to that value which requires Y 2 , = X ^tid 
Y, = |l,...,mj. Then the problem reduces to the one given 
in Chapter II, s^ = 0 for all i. 

As i ranges from one to m, let e^^ > e^> . . . > ep be an 
ordering of the ratios ~ 'ij i. . That is, let 

fn-i'iLi.D, . j I , 



= 



e 




max 

max 



f I 
UilvC 

( niiLiD 

L 




U t,' 



a rxi 

1 - 



A 

If < e^, the optimal values of , Sj^ , . . . , and s^^ are those 
which correspond to the least value of K from among the 
sequence of problems stated in the above two paragraphs. 



The first problem in the sequence is to 
i n i m 1 1 e 

Sl- Yi 

Subject to O +c»r cf. J 
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where Yj, and Y 2 , are the first Y^ and Yj, formed. 
For this problem 

K = EjQj 

f i n , niAiD, T \ . 

B.- ] > 

E, = * ,L •> 

^ ( l-T** 



icCi.4 



Lt-Yi 

.-On.-Dv ^ liCjllBl 

l«Y, i * 



Now K is convex in , and, hence, minimum K corresponds to 



Qi = 
Qi = 




otherwise 



e. such that KCe-,) - 



This procedure is followed p-1 times for the p-1 inter- 
vals for Qi , and for each iteration a value for minimum K in 
Q i is generated. The value for K where s^ = 0 for all i need 
not be considered since the last of the sequence of p-1 
problems examines K(ep) . From among these p-1 alternatives, 
pick the least value for minimum K. This will give a feasible 

^ /s 

value for Qt from which feasible s^ can be computed for all 
i e Y^ of that iteration. Of course, s^^ = 0 for all i £ Y^ 
of that iteration. 

CASE THREE : = 0 and > 0 . 

In this case equation (5-8) cannot be used to solve for 
Q since ^ | defined. Let > 0 be arbitrarily 

large but fixed. Then s^ is given by equation (5-7), and the 
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positive definiteness of K in Sj,,...,Sj^ implies that, for 
fixed , Sj , . . . ,Sjjj provide a global minimum. 

Now consider K as a function of Qj, at the point 

(s^ , . . . , • 



K 



= ^ fexpreS5.'or^C5-‘1)j + C 21 Dl + . 



(5-11) 



For Qj_ ^ e^, equation (5-11) is valid and minimum 
K = livYN K - C , 

G?i CO 

But suppose that Qj^< e . Then equation (5-11) is no 

longer valid, and an analysis similar to that performed in 

CASE TWO is required. The only modification is that E ^ 

Tw C D ■ 

reduces to E . The minimum K is picked from 



Y. 



5 D 



among the p-1 alternatives. If that value of K is less 
fhan C.LDt+'rit’iD:, then ,Sj^ , . . . ,s^^ are the corresponding 

A 

values for which that K applies. Otherwise, =oo and 

A 

s^ = <?0 for all i. This suggests it is best not to go into 
business . 



CASE FOUR ; < 0 and Bj^ = 0. 

In this case equation (5-8) cannot be used to compute 

A A 

since it implies that = 0 which in turn would imply 
s^ < 0 for all i. As before, minimum K for a given large 

A 

can be found when s^ is determined by equation (5-7). 

Then K = CLDl + X'TLv'Di + Q, (reC L \ • Since K is 

strictly increasing in Qi., K is minimized as Q i is minimized. 
But this argument is only valid for Qj, ^ . Hence one 

candidate for minimum K is 



I! 






(5-12) 



_ _ C— . ^ (ni-On.Dl 

K = C It>i + S, [loC L 2 n,i>\ 



For Qj_ < = 0 for some i; and the formula for minimum 

K changes. 

Since = 0 for some i, B increases; and, hence, 

B 2 > 0. Now CASE TWO applies. The minimum K under CASE 
TWO conditions is computed and compared to equation (5-12). 
If K in (5-12) is less than the CASE TWO solution, = e^_ 
and s^ is computed with equation (5-7). Otherwise, the 
CASE TWO solution applies. 



CASE FIVE 



Bj^ < 0 and Bj < 0 



In this case equation (5-8) cannot be used to compute 

f 

since \ is imaginary. As before, for a given 

large value of Qi the minimum value of K is given by: 

Ka ^ |exprcssic>n cZl>c + r^^'>.^^i+ Qi|roCZ^-~7~^'’| 

(5-13) 

Since expression (5-9) is negative, K is minimized for min- 
imum . But equation (5-13) is only valid for Qt ^ ej^ . 
Hence, a candidate for minimum K is the K given by equation 
(5-13) evaluated at = e^. 

F'or Qj_ < e^ construct and as was done in CASE TWO, 
Here, K is given by 

I f Aol>. . j «L^'_ T ^ ^ j-j.. 






_ . ^ f iLCn.Dll 

^ L^v3)c Q, jloC.L 2«;1>, ^ 2n;I>, J. 



(5-14) 



46 



Now, if 




remains negative, then the next candidate for minimum K is 
K given by expression (5-14) evaluated at e^. This analysis 
continues (evaluation at e^ would be next) as long as 
epxression (5-15) remains non-positive. If expression (5-15) 
turns positive at some iteration, then the procedures of 
CASE TWO apply for the remainder of the analysis. Then 
minimum K is picked from among the p alternatives, and 

^ \ A 

Qj_ ,s^ , . . . , and s^ are the values that pertain to the K thus 
picked . 

Although this search is sufficient to find minimum K 
it can be streamlined somewhat so that fewer alternatives 
need be considered for minimum K. As K is transformed from 
iteration to iteration and indices migrate from Y 2 , to , 

K is continuous in Q^. That means that K for two successive 
iterations when evaluated at the common boundary point is 
the same. While expression (5-9) remains negative we have 
seen that minimum K in is taken on for equal to the 
left hand or smaller boundary value. Thus, as K changes 
from iteration to iteration and as long as expression (5-9) 
remains negative, K continues to decrease with decreasing 
Qi. Accordingly, past values of K may be discarded as 
potential minimum K since they are greater than the present 
K. 

The first value of K that should be retained is the last 
value of K for which expression (5-9) is negative which 
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equals K evaluated at the right hand or upper boundary value 
for Qj^ for which iteration expression (5-9) is positive for 
the first time. 

CASE SIX : = 0 and < 0. 

This case is the same as CASE FIVE except that the term 
[loC. L I missing. If expression (5-15) 

turns positive, then this problem falls into CASE THREE 
except that = 21 appropriate 

W V ^ ^ I 

A 

^ \ \ 

values of Qj^ , s^ , . . . , and s^ are found as they were before. 

C. THE SEARCH FOR OPTIMAL n. , . . . ,n 

*■ m 

Finding optimal n^,...,nj^ involves many calculations. 

The first problem that must be addressed is, as in the 
no-stockouts-allowed case, showing that only finitely many 
points need to be considered. In general this cannot be 
done since, without some restriction on the parameters, as 
the values of n^^,...,n^^ range, the solutions for 

^ ^ A 

Q 4 ,s. , . . . , and s found in Section B cannot be controlled 
and the value of minimum K cannot be predicted. However, 
if the assumption that ZAiTiC ^ TCt^ D l for each i is 
made; the hyperbox G can be constructed. 

Consider the following two problems: 

(1) minimize K(q^,s^ , . . . ,s^^,n^ , . . . ,n^^) 
subject to 1 integer; 

Q 1 . ^ 0 ; and 

O < Qi -for all 
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I 



( 2 ) 



minimize KCQ^jS^,.. 
subj ect to ^ 1 

Q t ^ 0 j 



• 9 i 9***9 ^TT» ) 

m •*■ m 

integer; and 

and 

Qi for all i. 



The first problem is the minimization problem that CASES ONE 
through SIX address with n^,...,njj^ fixed. The second problem, 
given that 2/JiriC ^ has as a solution s^ and 

given by equations (5-7) and (5-8) for a given n^,...,nj^. 

The second problem is also a relaxation of the first problem. 
Let * (n^ , . . . ,n^^) be the optimal solution to the first 
problem for a given nj_,...,nj^. Let (nj^ , . . . ,n^^) be the 
corresponding solution for the second problem. Then 
Kl ^ ~ K;L (|^,1 • ^ •^m) , 

The following theorem shows that n* = (1,...,1). 

Theorem 5-2 : 2A*^.I{C for all i, 

then n* = (1 , . . . , 1) . 



Proof : The following argument shows that, if 

for all i, then whenever n^ ^ 2, regardless of the choice 

of n^ , 

K* K^*(l,...l) < K 2 ,*(nj ,...,nj^ (n^ , . . . ,nj 

where K* = (n^* , . . . ,n^^*) ; and, hence, n^* = 1 for each i. 

Substituting Qj_ given by equation (5-8) and s^ given by 
equation (5-7) into equation (5-4) for K yields equation 
(5-16), the solution to the second problem for given 

* • * »^m‘ 
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where 



Bj - 21 X *) 



BgCn) = 



B^(n) = 






Note that i£ n. ^ 2 for some i (or Bj_ < 0, the condition 



that makes this solution for K valid), then Bj, B^Cn), B^(n) , 
B^, B^ (n) , and Bg (n) are all positive. Then, for K given 
by equation (5-16), 



Consider now K(l, . . . ,1) . Since B^ > 0 and Bj^ = 0, 
CASE THREE applies. Recall that for CASE THREE 

Min KCQ.) ^ '^*'1 KCQ.)- 3- 

It follows that whenever 5 . 2 for some i 



This completes the proof. 

D. THE CONVEXITY OF K (Q^ , s , . . . , s^^^) 

Throughout this chapter the various efforts tihat have 
been made to solve for minimum K have disregarded the notion 
that K may not be convex in ,s^ , . . . , and s^^^ . In general. 



) > 63 . 




1) ^ Kj^*(n) < Ki*(n) . 



50 



it is not. One way to show convexity is to show that the 
Hessian of K, the matrix of second partial derivatives, is 
positive definite. The Hessian of K is 









3>^k 




30,dS| 


S0.3Sj. • • * 


^QidSrv% 


y>^K 




3^K 




35,‘^Q, 


ss,** 




e 







• • • 






3s ‘ 



Now, 



+ “L. Cy j j TCiSync 

1 ^ ^ ^ ^ n. 



. lAoD, 

= - T-. C s i TU . 

" n.Q.'*- ' 

■ O » 

Let (h..) be the Hessian of K. Since v _ 

13 oS;3Sj 



iD, 



=r O 



for i ^ j, (h^j) is sparse (if m = 9, then (h^j) has only 
28, out of 100, non-zero entries). However, the analysis of 
(hij) for a problem of arbitrary size is difficult. The 
case for m = 1 illustrates this. 
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<t?mp 



m- 



Suppose m = 1. Then has only four entries and is 

positive definite provided that hj^ > 0 and j > 0. 

For Q. > 0, h,| > 0. Now ^ 

. , I 1 . ^ fl.CAo ^ T 2 I|Cs,tCi 

h..Kxi-h,i - 0.-^1“^ xx>. T), 

If = 2, this reduces to 

* A.I.C-x,’- 

©.•^ I n. 

which is negative for sufficiently large. Hence, C^^j) 
is not positive definite for all feasible values of and 
s^; and K is not convex. 

Although the convexity of K is a sufficient condition 
to claim that a feasible stationary point of K is a global 
minimum, it is not a necessary condition. The analysis in 
Section B showed that if K was first minimized over 
Sj^ , . . . , and s (K is convex in s . . . , and s ) and if the 



resulting cost function was then minimized over 

(K(Q^ ,s^ , . . . , Sj^) is convex in Qi ) , a global minimum could 

be found. The following theorem formalizes that approach. 

Theorem 5-5 : Let 

KCQ..S'> 

and let ^ 

Kfo ) = vv»iio 1^ CQO • 

CLx 

Then K(Q^*,s *,..., s^*) = K(Qi) ; Qi* ^ Qi. ; and s.* = s.(Qi) 
for each i. 

Proof : Since KCQi* , . . . ,Sj^*) is the global minimum, 

K(Qj^*,s *,...,s *) < KCQi). 
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Now, 



K«0 S 
and, hence, 

K(Ql) ^ K(Q^*,s^S...,s/). 

Therefore , 

K(QJ = K(Qi*,s 

Finally, it follows that r Qt ; and si* = s^(Q^) for 

each i. This completes the proof. 

This method of finding the optimal solution is used 
throughout CASES ONE through SIX. It suggests a weaker 
sufficiency condition for (Q^* ,s^* , . . . ,s^*) than convexity 
of K. This condition is proposed in the following corollary 

A 

Corollary 5-4 : Let s^ be the value of s^ which minimizes 

A 

K(Qi , Sj^ , . . . , Sj^) for a given value of Qi. . Then KCQ^) = 
K(Qi,Sj^(Q^) ,. . . ,s^(Q^)) is convex. 

E. CLOSING REMARKS 

This chapter has offered the development of a backorders 
allowed extension of the model presented in Chapter II. 
Section B derived a method that can be used to find candi- 
date points for a local minimum to K. Section C showed that 
n.* = 1 provided that 2A.I.C ^ Tl’’. D. for each i. Section D 
showed the solution technique discussed in Section B indeed 
provided the global minimum. In all, the general backorders 
problem has been solved except for the problem of bounding 

TT^ 

the search in n when 2A^I^C ^ some i. 
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VI. CONCLUSION 



A. SUMMARY 

This paper has presented two deterministic multi- echelon 
inventory models. The assumptions of the first model led to 
the development of an extension of the EOQ (no-stockouts- 
allowed) model. Through the presentation of a handful of 
theorems, it was shown how the optimal policy could be found 
either by exhaustive search in the completely general case 
or by a rapid steepest descent search in the case where the 
ordering cost is the same for each lower echelon activity. 

From the author's point of view the three most inter- 
esting aspects of this paper are the assumptions for the 
steepest descent algorithm, the effect of 

optimal n in the no- stockouts - allowed model, and the general 
method for finding optimal K(|i) in the backorders model. 

The first point of interest is the way in which the 
assumption A^^ = A^ for each i affects the search for optimal 
n in the no-stockouts -allowed model. That this assumption 
streamlined the search is of some interest mathematically, 
but it is of additional interest from an economic point of 
view. In military applications of inventory theory it is 
usually assumed that the cost of ordering or requisitioning 
is constant for any one class of activities. For example, 
among Navy stock points inventory "pull" item orders have 
roughly the same cost. In the case of afloat activities 
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it was long assumed that the cost of submitting and pro- 
cessing any requisition was $7. The exact figures are not 
important here, but the principle lends additional justifi- 
cation for the assumption. 

The second point of interest is that the observation 
that implies that n^* = 1. This says that the 

cost to the system of holding goods at the upper echelon is 
at least as great as the cost of holding them at that par- 
ticular lower echelon. In this case, the solution shows 
that the optimization strives to minimize the holdings of 
the upper echelon because of the relatively high "opportunity 
cost" presented. This notion seems to support military 
decisions to provide central procurement activities with 
an unencumbered "open to buy," because the "opportunity 
costs" for alternative decisions are high. 

The third point of interest is strictly a mathematical 
one. Even though the cost function in the backorders - allowed 
model is not convex in all the continuous variables simul- 
taneously, theorem 5-3 and its corollary show a global 
optimum can be found by decomposing the cost function 
minimization problem into a sequence of two convex-program- 
ming problems. Thus, for this problem, a weaker sufficiency 
condition than the standard one of convexity has been found. 

B. RECOMMENDATIONS FOR FURTHER WORK 

Specifically, there are two points in this paper that 
should be re-examined and, perhaps, extended. The author 
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feels that an optimal search plan, such as the algorithm in 
Chapter III, can be developed for the general no- stockouts - 
allowed model. For the backorders allowed model it was 
shown that if ^ Tt^ for all i that n^* = 1 for 

all i. A bound and search scheme should be developed for 

the case where 2A.I.C< TL . D. for some i. 

11 11 

A general area suggested for further work is to extend 
the backorders- allowed model for the case in which the 
backorder-penalty cost may be time-weighted. Additionally, 
a "lost sales" model could be developed. 

C. CLOSING REMARKS 

In appraising the value of this new model it should be 
remembered that it represents an attempt in basically a new 
direction. Simple results usually forerun more complex ones 
which are often more directly useful. However, just as the 
standard EOQ model serves as a good first approximation to 
more complex deterministic or stochastic single activity 
models; so perhaps this model may be found similarly useful 
for multi-echelon models. 
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